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Abstract 

The dynamic susceptibility Xq{^) of the isotropic XY- model (s=l/2) on 
the alternating superlattice (closed chain) in a transverse field h is obtained 
exactly at arbitrary temperatures. It is determined fi'om the results obtained 
for the dynamic correlations < Sj,^{t)S[^{0) >, which have been calculated 
by introducing the generalized Jordan- Wigner transformation, by using Wick's 
theorem and by reducing the problem to a diagonalization of a finite matrix. 
The static properties are also reobtained within this new formalism and all 
exact results are determined for arbitrary temperatures. Explicit results are 
obtained numerically in the limit T = 0, where the critical behaviour occurs. 
A detailed analysis is presented for the behaviour of the static susceptibility 
Xq^(O), as a function of the transverse field h, and for the frequency dependency 
of the dynamic susceptibility Xgi^)- It is also shown, in this temperature 
limit, that within the magnetization plateaus which correspond to the different 
phases, even when the induced magnetization is not saturated, the effective 
dynamic correlation, < ^ ^jnW'^imi^) > ' time independent, which 

constitutes an unexpected result. 
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1. Introduction 

The homogeneous one-dimensional XY-model introduced by Lieb et al.[l] 
has been object of many studies. Its importance comes mainly from the fact that 
it is among the few many-body problems which can be solved exactly. Although 
the static and dynamic properties of the general homogenous model have been 
studied since its introduction (see e.g. Ref. [2] and references therein), the 
study of the inhomogcncous periodic model has boon restricted to the alternating 
chain [3-5] . Only recently have more general inhomogeneous models been object 
of study , namely, the alternating superlattice [6-8] and the inhomogeneous 
periodic model composed of inhomogeneous segments of finite length [9,10]. 
The results of these papers have been restricted to the static properties of the 
isotropic model and up to this date, to the best of our knowledge, no results 
have been obtained for the dynamic properties of these models. 

In this paper we study the dynamic properties of the isotropic model on the 
alternating one-dimensional superlattice. The dynamic properties are restricted 
to the longitudinal direction and they correspond to an extension of our previous 
results [6-8]. In section 2 we present the solution of the model and the basic 
results. In section 3 we obtain the dynamic correlations in the field direction 
and in section 4 the longitudinal dynamic susceptibility. Finally in section 5 we 
present the results and the main conclusions. 

2. The model: basic results 

We consider the isotropic XY model on the one-dimensional alternating 
superlattice, which we have been able to solve exactly [7]. The superlattice 
consists of N cells, composed of two subccUs, A and B, with ua and ub sites, 
respectively. The l-th unit cell is shown in Fig.l, and the distance s between 
two consecutive sites is taken equal to one. 

If we assume periodic boundary conditions for a chain with N cells, the 
Hamiltonian of the XY-model [7] can be written in the form 



H 



'UA — I 



1 ^ + - + - 

1=1 \ \_m=l m=l 



l.m ( ' 



(1) 



m=l 



m=l 



where I identifies the cell, S"^ = ± iS^, J is the exchange parameter between 
spins at the interfaces, Ja{Jb) the exchange parameter between spins within 
the subcell A{B), and hA{hB) is the transverse field within the subcell A{B). 
The spin operators can be expressed in terms of fermion operators using the 
generalized Jordan- Wigner transformation [6] 
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l — l riA 



l-l 



n=l r=l 



n=l r=l 



(2) 



I riA 
71=1 r=l 



r=l 



n=l r=l 



(3) 



and, by introducing this transformation, the Hamiltonian can be written in the 
form 



N (UA-I 



TLB -I 



1 = 1 m=l 



nA UB ^ 

+ E 2/*^ (4m«;,™ - l) + E 2/^3 - m + (4) 



where a's and h's are fermion operators, and given by 



AT riA 



N UB 



i=l r=l 



(5) 



is a boundary term which will be neglected. As it has been shown [11], this 
boundary term, in the thermodynamic limit, does not affect the excitation spec- 
trum, the static properties of the system, nor the dynamic correlation function 
in the field direction. Introducing the Fourier transforms [8], 



BQ,m = E {-iQdl) b. 



,m 5 



(6) 
(7) 
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where Q = 2-1171/ Nd, n = 1,2, N, and d = UA + ns is the size of the cell, the 
Hamiltonian can be written as 



Q 

where Hq is given by 

nA-l / J \ riB-l / J \ 

E (^^k™^Q,-+i+/^-cj - E [-fBl^BQ,^+, + h.c^ - 
^ K-^^>^ + <n.^<^^e-^' + + ivi!M^^^±!^. (9) 

We can also write H in the form 

H = Y,V^T{Q)Vq, (10) 

Q 

where 



^^=(4,i'-'<nV<i'-'^Q,.B) (11) 

and T((3) is a matrix of dimension x which represents the quadratic Hamilto- 
nian Hq. Since [iJg^ J^q/] = 0, the diagonalization of the Hamiltonian is reduced 
to the diagonalization of Hq, which can be written as a free fermion system 



-'^Q = E ^Qk^Q,k^Q,k^ (12) 

k 

where Eq). are the diagonal elements of UqTqUq, and Uq is the unitary 
transformation which diagonalizes iJg, which is determined numerically. The 



N 



V 2 



(8) 
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energy spectrum presents d branches and, by using the unitary transformation 
Uq, we can express the fermion operators a's and b's in terms of the ^'s, as 



k 

'^Q,nA+m,k^Q,k^ 

k 

-j= exp {iQdl) UQ^m,kiQ,k, 

Q,k 

bl,m = -y^Ye^P{iQ<ll)UQ^nA+m,k^Q. 
Q,k 



O-Lm 



k' 



(13) 
(14) 

(15) 

(16) 



where UQ^kj are defined as 



UQ,k,l 

UQ,k,nA 
UQ,k,nA+l 

1^Q,k,nA+nB 



Eqk 



UQ,k,l 

'^Q,k,nA 
UQ,k,nA+l 

UQ^k,nA+nB 



(17) 



The general solution of this equation is determined numerically, although 
analytical solutions can be found for some special cases. For instance for ha = 
ub = ^, JiA = JiB = h, we can express explicitly the solution in the form 



u;Q = -h± -^^c±v/^(Q^ , (18) 



where 



c = 



2 



{A + Jl) 



g{Qd) 



\j^JaJb cos m + ^ ( Jl + Jl) + ^ {Jl 



JV 



(19) 
• (20) 
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As it has already been shown [7], the effect of a homogeneous field, Ha = 
/iB = ft, is to shift the zero field spectrum. For the special case shown above 
this can also be shown directly from eq.(18) and, consequently, the existence of 
a mode of zero energy will depend on the strength of the field. 

In passing, we would like to note that the spectrum can also be calculated 
exactly by using the position space renormalization group approach [12], and 
approximately by using a transfer matrix method [13]. Although the latter is an 
approximate method, we have shown that it reproduces numerically the exact 
result. 

Within the formalism introduced above we can obtain easily all the thermo- 
dynamic properties of the system. Then, we can express the internal energy U 
in the form 



U = J2^Qk i^Qk) , with {HQk) = ^ l^pj^^^ 

Q,k 



(21) 



where /? = l/ksT, Ub is the Boltzmann constant and T, the absolute temper- 
ature, and from this expression we obtain immediately the specific heat C 



(22) 



Q,k 



The induced magnetization per site and cell, which is defined as [7] , 



1 



Li=l \rn=l 



m=l 



(23) 



can also be written in the form 



N {nA + ub) 



Qk 



1 

2' 



(24) 



and from this expression, by making the identification Ha = = h, we can 
obtain the isothermal susceptibility, Xt — M i'^cei) > which is given by 
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^T= j,J^,^^ I.M'-MPE^,). (25) 



3. The Dynamic Two-Spin Correlation Function in 
the Field Direction 

The dynamic two-spin correlation function, {Sj^^{t)Sf_^j^^^{0)) , can be eas- 
ily obtained by writing explicitly the time-evolution of the operator Sj^^{t). 
This is obtained by using eqs. (2) and (15), which gives 



S^it) = ^ Yl [e^PH{Q-Q')dj)exp{i{EQ,-EQ,kOt) 

Q,k,Q',k' 

'^*Q,k,m'^Q' 'k' ,m^Q,k^Q' ,k'] " \, (26) 



and from this we can write 



SfUt)StlR,n{^)) =< f ^ E (-^W - Q')dj) X 

y Q,k,Q',k' 

X exp{i{EQk - EQ'k')t)u*Q^k,mUQ',k',m^Q^k^Q',k' - ^ 

^[j^ E exp(-i(Q-g')(iOwQ,fe,«WQ'.fc',™CQ,fe?g',fc' - ^ > • (27) 

\ Q,k,Q',k' J 



By using Wick's theorem we obtain 
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(0) 



^ exp {-iQRdl) exp{iEQkt)u*Q^k,mUQ„k„m {nqk) j x 
^'^exp{iQRdl)exp{-iEQkt)u*Q,.,^UQ,,k„n (1 - («Qfe)) J + 



Q,k 



\ Qk J \ Qk 



UQ,k,nU*Q^k,n (nQk) + 



'^Q>k,mU*Q^k,m {nqk) '^Q'k,nU*Q,k,n {^Qk) + \ , (28) 



Qk 



2N 



Qk 



which can be determined numericaUy. 

We can also obtain the dynamic correlation between average cell spin oper- 
ators, rf , defined as 



{ha + Ub) 



(29) 



m=l 



which satisfy the equality (rf) = (S^^i) . Therefore, by using eqs. (28) and 
(29), we can determine the dynamic correlation function ("rf (f)rjl|_^(0)) by ex- 
pressing it in terms of dynamic correlations between site spins and it can be 
writen in the form 



(rr(i)rf+«(0)) = 



NinA 



Q,fe 
2 



+ 



{UA + ub) 



exp {-iQRdl) exp{iEQkt) (nqk) j x 



.Q,k 



^ exp (iQRdl) exp{-iEQkt) (1 - {nqk}) J • (30) 



4. The Longitudinal Dynamic Susceptibility 
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The dynamic susceptibility Xq{^) can be determined by introducing the 
spatial and temporal Fourier transforms 



(T^(i)riQ(O)) =^exp(-iQi?)(Tf(i)rf+^(0)), (31) 



and 



1 r°° 

(t^qT^q)^ = - y ^ expM) (r^(i)rlQ(O)) dt, (32) 



and by using eq. (30) we can write 



+ 



N{nA + nB) 

X (1 - ("Qi-Qfe)) • 



(33) 



[14] 



From this result we can obtain the dynamic susceptibility from the expression 



Xq (w) = -27r « t'q; t^q » 



-f 

J — C 



(1 - exp(/3a;)) {tqt'Lq) , duj' 



(34) 



which can be written as 



N {ha + nsf ^ + ^Qk - Eq^-Qk 



E 



(35) 



In the limit w ^ 0, Q — > 0, we obtain the isothermal susceptibility as in the 
uniform model [15]. 
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The real and imaginary parts of Xq {^) can be determined from the previous 
expression by considering Xq ('^ " i^) in the hmit e — > 0, and from this we 
obtain expUcitly the results 



^Xq{<^) = 



N {ua +nB) 



1 



Qi,k 



(36) 



and 



TT 



N{nA + nB) 



Ql,k 



(37) 



5. Results and Conclusions 

The results for J ~ I, Ja — 2, Jb = fiA = fiB = h and ha = ub = 2, 
at T = are shown in Figs. 2, 4, 5, 8, 10 and 12. For the same interaction 
parameters and ua = 2,nB = 3, the results are shown in Figs. 3, 6, 7, 9, 11 
and 13 . 

In Figs. 2 and 3 we present the magnetization (rf) and the isothermal 
susceptibility Xt functions of the field. In these cases, as expected, we have 
four and five critical fields respectively, associated to the quantum transitions 
induced by the field and the isothermal susceptibility diverges at these points. 
In the first case {n = ua + ns, even) there is a plateau at zero magnetization, 
whereas in the second one (n = UA+riB, odd) the magnetization plateaus are at 
non-zero values. These results suggest that these plateaus, which are associated 
to the gaps in the energy spectrum [7], can also be related to the cluster spin 
states in a representation where the effective spin is diagonal [10]. 

In Figs. 4-7 we present the real and imaginary parts of the dynamic corre- 
lation function (Tf (i)Tf^^(O)) , as a function of time t, for R = 1. For n = 4, 
Figs. 4 and 5, the results are presented for three different values of the field, 
namely, h = hjc and h = hjc -I- £ (e = 0.001, 0.01), where j = 1, 3 respectively 
for each figure. For n = 5, Figs. 6 and 7, the results are also presented for three 
different values of the field, namely, h = hjc and h = hjc — e {e = 0.001, 0.01), 
where j = 1,3 respectively for each figure. In all cases, the imaginary part of 
the correlation goes to zero as we approach the magnetization plateaus, and the 
real part tends to the square of the magnetization (rf) . These results show that 

within the plateaus, although (Sj!^ i^)^j+R ni^)/ '"^^^ '^cit constant, the field 

does not induce any quantum fluctuation on the correlation (^f (i)Tj^_|_^(0)) . 
This means that the system, although not saturated, is in a frozen effective spin 
state. 
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In Figs. 8 and 9 wc show the static correlation function (r| as a 

function of R, the distance between cells, at T = 0, and three different values of 
the field. The correlation shows an oscillatory behaviour which is present for any 
value of the field outside the plateaus, with increasing period as the critical point 
is approached. On the other hand, at the critical field there is no oscillatory 
behaviour, which means that the period of the oscillation tends to infinity, and 
this result is still valid for any value of the field within the plateaus. This is 
consistent with the scaling form and the analytical continuation proposed for 
the two-spin correlation function (Sj (0)5;^(0)) in the homogeneous XY-model 
[16], where the correlation length is associated to the oscillation period. 

We have also verified that for h > hie {n = 4, even) and h > h^c {n = 
5, odd), since the system is saturated, the dynamic correlation function < 
^fn^^^ ^^)^im^^^ docs not depend on time nor field, while in the inter- 
mediate plateaus it presents a time dependence which is independent of the 
value of the field. 

The static susceptibility Xq^(O), at T = 0, is shown in Figs. 10 and 11 for 
various wave- vectors and the two different cell sizes. As already pointed out, 
the susceptibility Xo^ (0) is identical to the isothermal one Xt" ^'^'i consequently 
diverges at the critical points. Also, independently of the wave- vector , the static 
susceptibility goes to zero within the magnetization plateaus since in this region 
the spacial modulation of the magnetic field does not induce any fluctuation in 
the magnetization. The singularities present for the different non-zero wave- 
vectors are related to the oscillations of the correlation (^f (i)Tf^jj(O)) , whereas 
the ones for Q = are related to the critical points. As pointed out by Jullien 
and Pfeuty [17] for the homogeneous model and Lima and Gongalves for the 
model on the supperlattice [18], the non-critical singularities can also be asso- 
ciated to the unstable critical points. which appear in the study of the system 
within the real space renormalization approach. 

Finally, in Figs. 12 and 13 we present the real and imaginary parts of the 
dynamic susceptibility, xl^{uj), at T = 0, for values of the field close to critical 
ones and outside the plauteaus, and the same parameters of the previous figures. 
As expected, for a given wave- vector, the bandwidth of the imaginary part of 
the response depends on the field strength and on the size of the cell and there 
arc no infinite singularities. Also as expected, to the finite discontinuities at the 
band edges of the imaginary part correspond divergences in the real part. It 
should be noted that for values of the field within the plateaus the response 
goes to zero for any non-zero wave-vector. 
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Figure Captions 



Fig.l-Unit cell of the alternating superlattice. 

Fig. 2- Average magnetization per unit cell (a) and the susceptibility in the 
field direction (b) as functions of the field, for ua — tib = 2, 
J = 1, Ja — 2, Jb — 3, and Ha — fiB — h, a,t T = 0. The critical 
fields are hd = 0.691, hc2 = 1.096, and /ics = 1.5960. 



Fig. 3- Average magnetization per unit cell (a) and the susceptibility in the 
field direction (b) as functions of the field, for ua = 2, = 3, 
J = 1, Ja = 2, Jb = 3, and Ha = fiB = h, a.t T = 0. The critical 
fields are hd = 0.207, hc2 = 0.935, /ics = 1-207, hd = 2.161 
and hc5 = 2.226. 



Fig. 4- The real (a) and imaginary (b) parts of the correlation function 

(T|(t)T|_|_j(0)) , as functions of the time, at T = 0, for ua = ns = 2, 
J = 1, Ja = "2, Jb = 3, and h = hic (plateau), h = hic + £, 
with £ = 0.01 and 0.001. 



Fig. 5- The same as in Fig. 4 for /i = h^c (plateaux), h = h^c + £, with 
£ = 0.01 and 0.001. 



Fig. 6- The real (a) and imaginary (b) parts of the correlation function 

{'i~j{t)'i~j+i{0)) 7 as functions of the time, at T = 0, for ua = 2, ns = 3, 
J = 1, Ja = "2, Jb = 3, and h = hic (plateau), h = hic — £> 
with £ = 0.01 and 0.001. 



Fig. 7- The same as in Fig. 6) for h = h^c (plateaux), h = hzc — £, with 
£ = 0.01 and 0.001. 



Fig. 8- The static correlation function (r| Tj+fl) as a function of R (distance 
between cells), at T = 0, for = = 2, J = 1, = 2, Jb = 3, 
and for values of the field near and at the critical field. 
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Fig. 9- The static correlation function (t| t|_|_^) as a function of R (distance 
between cells), at T = 0, for ua = 2, ub = 3, J = 1, Ja = 2, = 3, 
and for values of the field near and at the critical field. 



Fig.lO- Static susceptibility in the field direction, Xq^(O), at T = 0, as a 

function of the field for ha = ns = "2, J = 1, Ja = Jb = 3, and 
Ha = hs = h, and different values of q. 



Fig.ll- Static susceptibility in the field direction, Xq^(O), at T = 0, as a 

function of the field for ua = 2, ns = 3, J = 1, Ja = 2, Js = 3, and 
Ha = fiB = h, and different values of q. 



Fig.l2- The real (a) and imaginary parts of the dynamic susceptibility in 
the field direction, Xq^{^)> at T = 0, as a function of frequency for, 
nA = nB = 2,J =1, Ja = 2, Jb = 3, and h = 0.693. 



Fig.13- The real (a) and imaginary parts of the dynamic susceptibility in 
the field direction, Xq^(a;), at T = 0, as a function of frequency for, 
riA = 2,ns = 3, J = 1, J a = 2, Jb = 3, and h = 0.2. 



14 



This figure "FIGl.JPG" is available in "JPG" format from: 



http://arXiv.org/ps/cond-mat/01 1 1372vl 



This figure "FIG2.JPG" is available in "JPG" format from: 



http://arXiv.org/ps/cond-mat/01 1 1372vl 



This figure "FIG3.JPG" is available in "JPG" format from: 



http://arXiv.org/ps/cond-mat/01 1 1372vl 



This figure "FIG4.JPG" is available in "JPG" format from: 



http://arXiv.org/ps/cond-mat/01 1 1372vl 



This figure "FIG5.JPG" is available in "JPG" format from: 



http://arXiv.org/ps/cond-mat/01 1 1372vl 



This figure "FIG6.JPG" is available in "JPG" format from: 



http://arXiv.org/ps/cond-mat/01 1 1372vl 



This figure "FIG7.JPG" is available in "JPG" format from: 



http://arXiv.org/ps/cond-mat/01 1 1372vl 



This figure "FIG8.JPG" is available in "JPG" format from: 



http://arXiv.org/ps/cond-mat/01 1 1372vl 



This figure "FIG9.JPG" is available in "JPG" format from: 



http://arXiv.org/ps/cond-mat/01 1 1372vl 



This figure "FIG10.JPG" is available in "JPG" format from: 



http://arXiv.org/ps/cond-mat/01 1 1372vl 



This figure "FIG11.JPG" is available in "JPG" format from: 



http://arXiv.org/ps/cond-mat/01 1 1372vl 



This figure "FIG12.JPG" is available in "JPG" format from: 



http://arXiv.org/ps/cond-mat/01 1 1372vl 



This figure "FIG13JPG" is available in "JPG" format from: 



http://arXiv.org/ps/cond-mat/01 1 1372vl 



